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*> Abstract 

<N 

We present a self-gravitating dyon solution of the Einstein- Yang- 

Mills-Higgs equations of motion in asymptotically AdS space. The 
back reaction of gauge and Higgs fields on the space-time geometry 
^) ■ leads to the metric of an asymptotically AdS black hole. Using the 

^D . gauge/gravity correspondence we analyze relevant properties of the 

finite temperature quantum field theory defined on the boundary. In 
particular we identify an order operator, characterize a phase transi- 
S^ . tion of the dual theory on the border and also compute the expectation 

H ' value of the finite temperature Wilson loop. 

1 Introduction 

The AdS/CFT correspondence P-[3] originally proposed to connect gauge 
and string theories has recently became a powerful tool to study strong 
coupling physics in systems of interest in quantum field theory and con- 
densed matter physics from the analysis of the low energy approximation of 
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string theories, namely gravitational theories [I]- [7] (for references on recent 
progress see for example [8]). 

In a previous work [9] we have used such gauge-gravity duality to relate a 
d — 3+1 Yang- Mills Higgs model in a Schwarzschild-Anti-de Sitter black hole 
background with a d = 2 + 1 quantum field theory defined on the boundary 
In this way we were able to study the strong coupling regime of the theory 
on the boundary from the classical dyon solution found in the bulk, showing 
that the 2 + 1 system undergoes a second order phase transition and exhibits 
the condensation of a composite charge operator. 

It is the purpose of the present work to extend our study to the case in 
which there is back reaction on the space-time geometry so that the complete 
Einstein- Yang-Mills-Higgs (EYMH) system of equations of motion has to be 
solved, with the condition that the metric solution corresponds to a black 
hole in asymptotically AdS space. 

Self-gravitating dyon solutions of the EYMH equations in asymptotica- 
lly AdS space, regular both at the origin and at infinity, were constructed 
in [TD]-|TT] . It was found that, as the value of the dimensionless parameter 
GnHq grows, the metric function develops a minimum which approaches to 
zero at a critical value (Gn is the Newton constant and Hq the v.e.v. of the 
Higgs scalar). Above the critical value a singularity at the origin develops 
indicating the possibility of a threshold horizon. In fact, black hole solu- 
tions that are finite at spatial infinity and at one horizon have been shown 
to exist in systems like the one we are interested for the case of asymptoti- 
cally flat space [I2]-J15j. One of the purposes of the present work is precisely 
to construct such black hole solutions in asymptotically AdS^ space, having 
finite mass and non-zero horizon radius. Since the black hole naturally in- 
troduces a temperature we shall be able to analyze relevant properties of the 
finite temperature quantum field theory defined on the boundary using the 
gauge/gravity correspondence. 

It is interesting to notice that, in the context of the AdS/CFT corre- 
spondence, gauge fields in the AdSd bulk can induce a field theory on the 
boundary containing dynamical gauge fields for low dimensions (d < 4) |16j - 
[2"U] . It then becomes of interest to compute the expectation value of the finite 
temperature Wilson loop of the theory at the boundary which, according to 
the gauge/gravity correspondence [21]-[25], is related to the Nambu-Goto 
action in the bulk associated to the dyon solution. 

The plan of the paper is the following. In section 2 we present the model 



in the bulk M and propose a spherically symmetric ansatz that reduces the 
EYMH equations of motion to a a nonlinear system of coupled ordinary dif- 
ferential equations. Such ansatz corresponds to a <9M = S 2 x S 1 boundary 
(with the time direction compactified to S 1 ). We also consider the appropri- 
ate change of variables leading to <9M = R 2 x S 1 boundary which could be of 
interest in view of condensed matter applications. In section 3 we discuss the 
appropriate conditions on the fields and metric in order to have a solution 
giving a black hole in asymptotically AdS± space and the gauge and Higgs 
fields of a dyon. We construct numerically such a solution and discuss its 
properties. We then apply the holographic correspondence to analyze the 
behavior of the system defined on the border as the temperature changes, 
identifying an order parameter and characterizing the phase transition that 
the 2 + 1 system undergoes. Section 4 is devoted to the calculation of the 
Wilson loop and binding energy for the theory in the border. Finally, we 
present a summary and discussion of our results in section 5. 

2 The model 

As in ref.[9], we consider a gravity- Yang-Mills-Higgs system with gauge group 
SU(2) and the scalar field in the adjoint representation, in a 4 dimensional 
space-time with Minkowski signature (— , +, +, +). The action takes the form 



S = S G + S YM + S H = J d 4 x y/\g\ (L G + L YM + L h ) (1) 



with 



■>G 



(R-2A) (2) 
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Lym = -^ f; v f^ (3) 

L H = --D^H a D^H a - ~ (H a H a - H 2 ) 2 (4) 

Here k 2 = 8 n G^ with G^ the Newton constant; e the gauge coupling and 
A the cosmological constant (we take A < which corresponds for our con- 
ventions and in the absence of matter, to anti-de Sitter space). For simplicity 
we will consider the BPS limit A/e 2 = 0. The field strength F^ v (a = 1, 2, 3) 



and the covariant derivative D^ acting on the Higgs triplet H a are defined 
as 

F; v = d^Al - d v Al + e ahc A\A c v , D, L E a = d u H a + e ahc A\R c (5) 

The most general static spherically symmetric form for the metric in 3 
spatial dimensions together with the t'Hooft-Polyakov- Julia-Zee ansatz for 
the gauge and Higgs fields in the usual vector notation for internal degrees 
of freedom reads [TU]-[TT] 

ds 2 = -fi(x) A(x) 2 dt 2 + /i(x)" 1 dr 2 + r 2 dfl 2 

A = dt e H J(x) e r — dO (1 — K(x)) e v + dip (1 — K{x)) sin 9 ee 

H = H H(x)e r (6) 

where H sets the mass scale ([H ] = m 1 ) and we have introduced the di- 
mensionless radial coordinate x = e H r. 

From S 2 to R 2 

As stated in the introduction, we are interested in the physical system de- 
fined on the boundary <9M using the dual classical description of the gravity 
system governed by the action ([1]) in the 4 dimensional manifold M which is, 
asymptotically, anti de Sitter space. This means that after compactifying in 
the time direction, one has <9M = S 2 x S 1 . Now, following the approach in 
[23] one can trade such boundary with <9M = R 2 x S 1 by introducing a dimen- 
sionless parameter R through a change coordinates (£, x, 9, if) — > (r, y, x l , x 2 ) 
(or equivalently (t, x, 9, cp) — > (r, y, p, <p)), and then taking the R — > oo limit. 
As we did in ref . [9] , we consider the change of variables 

r = Rt 

y = {jo/R) x = r/(LR) 

x 1 = 2 RL tan(#/2) cos if = p cos</? 

x 2 = 2 RL tan(6 l /2) sin(/? = p simp (7) 

together with the field redefinitions 

A( y ) = A(x) , f(y) = R- 2 fa) , 

H(y) = H(x) , K(y) = R" 1 K(x) , J(y) = i?" 1 J(x) (8) 



After this change, the ansatz (J5J) becomes 
ds * = _/( y ) A(y) 2 dr 2 + L 2 ^ ' - 2 ° 



A = dr eH J(y) 
+ 



/GO 

ALRp 
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4(LiJ) : 



4 (Li?) 2 - p 1 



p 2 + 4(Li?) 2 " p ' p 2 + 4(Li?) 2 C3 



4Lfl " '4-*m 



p 2 + 4(Li?) 2 V# 



4(Li?) 2 -p 2 



4LRp 



2 ^P 



if = H H{y) 



4LRp 



p 2 + 4 (Li?) 



p 2 + 4 (Li?) 

4 (LR) 2 - p 2 



p 2 + 4 (LRf 



2 "p 



p 2 + 4 (Li2) : 



e 3 p d<p - e«o dp 



(9) 



and the equations of motion for the Einstein- Yang Mills-Higgs system take 
the form 



y f(y) 



1 
R 2 



3y 2 -K 2 H 2 \f( y )V 1 + V 2 



y 2 J'(y) 2 , J(y) 2 K(yf 



2 A(y) 2 f(y)A( y y 



vi+ J K K ^ y }l \ 



f(y) 2 A(y) 2 
f(y)A(y)K'(y)\ -= A(y) K(y) 



'K(y) 2 -± 



H{yf 



7o' 



Ay) 2 

7o 2 f(y) A(y) 2 



y 2 f(y) A(y) H\y) = 2A(y) H(y) K(yf 



f(y)A(y) \ y -^-\ =2J(y)K(y) 2 



where 



(10) 



Vi = 7 o 2 ^'(Z/) 2 + tV ^'(?/) 2 » ^ = 7o : 



Here we have defined the dimensionless parameter 
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In the R — > oo limit ansatz (JHJ) becomes 

ds 2 = -f{y) A{yf dr 2 + L 2 J^- + y 2 (dx l2 + dx 2 

f(y) v 

K(y) 

A = dr e H J(y) e 3 H — (eidx 1 — e 2 dx 2 ) 

±j 

H = H H(y) e 3 (12) 

Then system (fTUI) reduces to 

(y Kv)) = 3y 2 - k 2 H 2 ( f(y) v x + v 2 
yA'(y)^ K 2 H 2 [v l + ^M)A(y) 



y 2 J\yf , J(y) 2 tf(y) 



2 A(y)* f(y)A(y) 



2 



/(y) 2 A(y)2 



(/(y) i(y) A%)) = i(y) £(y) -^- + — ff (y) : 



2 



r lo z lo z f{y)A{yf 

y 2 f(y) A(y) H'{y) ) = 2A(y) H(y) K(y) 2 



f(y)A(y) |^^| =2J(!/)^ 



where now 



(13) 



* = Vf£ (14) 

3 The black hole-dyon solution 

In reference [H] we constructed self-gravitating dyon solution in asymptot- 
ically AdS± space which were regular both at the origin and infinity. An 
important feature noted in that work was that as the dimensionless parame- 
ter kHq grows, the metric function / develops a minimum which, at a critical 
value (kH ) c , corresponds to a zero of /. For values of nH > (kH ) c dyon 
solution which are regular at the origin ceases to exist; however, as explained 
in [12]- [13] one can hope that solutions that are finite at spacial infinite and 
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at one horizon can be found and this will be precisely the issue we shall an- 
alyze in this section. We shall construct here such kind of black holes dyon 
solutions for the Yang-Mills-Higgs system in asymptotically AdS^ showing 
that they have finite mass, non-zero horizon radius and, consequently, an 
associated Hawking temperature so that one can study finite temperature 
effects in the quantum field theory defined on the boundary. 

3.1 Boundary conditions and properties of the solution 

We shall look for a solution to ( fit)]) with a horizon at nonzero radius y = yh 
that we can chose, without loss of generality, to be yh = 1- At the horizon 
we impose the conditions: 

/(l) = and A(y),K(y),H(y)J(y)/(y-l) regular at y = 1 (15) 

so that 



f(y) = h(y - I) + 0[(y - I) 2 ] 




A{y) = oo + oi(j/-l) + 0[(y - 


l) 2 ] 


H(y) = h + h 1 (y-l) + O[(y- 


-I) 2 


K(y) = k + k 1 (y-l) + O[(y- 


-i) 2 : 



J(y)=j 1 (y-l) + 0[(y-l) 2 } (16) 

Concerning y — > oo, the asymptotic behavior should be 

f(y) = y* + F + ^ + --- 

y 

A(y) = l + --- 

\6(y)\ = H Q -^ + ^ + ^ + ..., ,61 (17) 

*M = j^ + 5^.+"- as) 

J(y) = Jo + - + • • • (19) 

Such behavior is consistent with the K equation of motion whenever the 
following condition holds 

— = v(v + l) (20) 

7o 



Equation (12TJJ) has two solutions, 



u ± = ~±^y/l + Am w ^ (21) 

where m w = eH is the gauge boson mass. Only the u + root gives an 
acceptable asymptotic behavior for K(x) so that one ends with 

K(x) = — = -r t , as x ->■ oo (22) 

As in [9] we define the field strength associated with the surviving U(l) 
symmetry in the form 

^ (1) - ^ (23) 

with the magnetic and electric fields, and the corresponding charges, given 
by 



cr i J K 



1 s ijk 

Q m = J d?x^vf ] B* = ~ (25) 

Q e = J (Pxy/^V^tf (26) 

Note that the spherically symmetric solution has a quantized n — — 1 mag- 
netic charge (in units of Air/e). 

Concerning the black hole temperature, it is given by 

r=^ z A(y h )f(y h ). (27) 

Now, taking into account that we have y^ — l and defining the dimensionless 
quantity 

T = 7H a (28) 

one has 

T = JWi (29) 

Air 

where a and fi are defined in ffTBT) . Since a and /i depend both on k and 

7o we have that T = T(k, 70). Note that T(0, 70) = (3/47r)7 . 
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3.2 Numerical analysis of the solution 

We have solved system ffTUl) with boundary conditions (fT5l) - (fT9l using the 
relaxation method. In this method the differential equations are replaced 
by a system of finite-difference equations on a one-dimensional grid with N 
points (typically iV=500 in our case). Then, the solution is determined by 
starting from an initial guess and improving it iteratively (see ref.[2S] for 
details). 

The system (flOl) is composed by two first-order differential equations (for 
the metric functions f(y) and A(y)) and three second-order ones (for the mat- 
ter fields K(y), J{y), and H(y)), so to implement the integration algorithm 
we need eight boundary conditions. At the horizon, y^ = 1, we choose 

/(I) = 0, 
J(1) = 0, 

/'(i) i(i) k\i) = i(i) k(i) (/?(i) 2 - ± + ij # (i) 2 

/ > (l)l(l)fl'(l) = 2i(l)H(l)A'(l) a (30) 

The last two conditions are derived by simply evaluating the differential 
equations for K and H at the horizon. 

For the boundary corresponding to y — > oo we choose for the numerical 
calculation a finite size L = 10, and check that the solutions are stable under 
variations of L. We impose at this boundary the conditions 

A = l, 
H=l, 

K oc y^ 1 , 

J=Jo (31) 

Following this procedure, we have searched for regular solutions for dif- 
ferent values of parameters k and 70. The corresponding solution depends 
on just one boundary condition, namely Jo- 

As is standard in the relaxation method, its effectiveness is primarily 
determined by the election of the initial trial functions, especially in this 
system that comprise eight first-order differential equations. We started by 
solving the system in the region of small 70 and k where the algorithm is 
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Figure 1: The three regions in parameter space. Only in region I the non 
trivial self-gravitating dyon solution exists. Region II corresponds to the 
"trivial" solution and no solution exists in region III whose frontier indicates 
the appearance of a second horizon. 

very stable and the election of trial functions is simple. Once we have found 
solutions in this region, we select them as trial functions for the equations 
with larger constants. As we will discuss below, for certain critical values of 
the parameters there are no non-trivial solutions to the differential system. 
Near these regions the algorithm becomes unstable and more mesh points and 
a careful election of the trial functions is needed in order to obtain solutions 
with the desired accuracy. 

Our analysis allows to identify, three regions in parameter space (k, 7 ) 
with a completely different behavior, as depicted in figure 1. Region I is the 
one where the gauge and Higgs field solution correspond to a dyon solution, 
i.e. a soliton with both electric and magnetic charge. Concerning the metric, 
it corresponds to an asymptotically AdS dyonic black hole. We shall discuss 
in detail this solution at the end of this section. 

In region II, on the right, only a "trivial" solution exists for which the 
gauge and Higgs fields take the form 



H{y) = l, K{y) = 0, J = J (1-- 

V 



(32) 



Concerning the metric functions, they correspond to an AdS dyonic black 
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hole, the result of the back reaction on the geometry 

~, , ~, . 1 2 M M-l- 1/i? 2 , , 

A(y) = l, /(v) = ^ + ^-y + j^— (33) 



Here M 



M=i+^+^ l (j 2 o+ : i) 



R 2 2 \ u R 4 



The reason for the existence of region II is well known: a too large Newton 
constant destabilizes the dyon solution (its mass grows and its radius de- 
creases). This effect is enhanced in the presence of a cosmological constant 
so that both mechanisms lead to the occurrence of a critical line after which 
the dyon becomes gravitationally unstable and only the trivial solution exists 
(see [I] and references therein for a more detailed discussion). 

Finally there existe a region III on the top of the figure, which is bounded 
from below by a critical line indicating the appearance of a second horizon. 
No solution, besides the trivial one, exists above this line. For a more detailed 
discussion on this issues see [T2]-[T3]). 

We have analyzed the solutions both for boundaries <9M = S 2 x S 1 (fixed 
R) and for <9M = R 2 x S 1 (limi? — > oo). Qualitatively the solutions are very 
similar, thus for definiteness, we shall present in detail the latter case. 

Figures 2 shows a representative dyon solution arising in region I. One can 
see that the Higgs fields rapidly attains its symmetry breaking constant value 
while the magnetic and electric fields (associated to K and J respectively) 
concentrate in a spherical shell starting at the horizon. 

Concerning the metric function /, one can see that as one increases k an 
outer minimum starts to develop (see figure 3) and moves downward, finally 
reaching zero and becoming an extremal horizon at some critical value n cr a 
phenomenon already found for gravitational monopoles in the case A = |12j - 
[T3] . This defines the critical line between regions I and III. The coefficient 
/i, the slope of / at the horizon has only a small variation with k, however 
the coefficient a^ that represents the function A at the horizon, goes from 1, 
in the absence of back reaction, to at k = k c . 

We conclude by noting that for small k the results described above for 
the solution in region II are qualitatively similar to those obtained in [9] the 
case in which back reaction is ignored and a Schwarzchild-AdS4 black hole is 
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Figure 2: Solution for the gauge field functions (K, J and scalar field H) for 
7o = 0.1, kHq = 0.5 and Jo = 5.5. The solution exists starting at the horizon 
y = Vh = 1. 

taken as a background. This can be explained if one assumes the existence 
of a weak gravity limit so that the k dependence is smooth and the solutions 
change continuously. However, even in this situation similarities are lost 
for certain ranges of parameters (k, 70), as the phase diagram presented in 
Fig.l exhibits a critical line that signals the existence of a second horizon, 
above which there are no non-trivial solutions. One can associate a Hawking 
temperature to this second horizon but since the bulk to be considered within 
the gauge/gravity duality would correspond just to the exterior of this second 
horizon, where only a trivial solution exists, no relevant physics could arise 
in this case. 



3.3 Holographic correspondence at finite temperature 



As explained in section 3.1, the Hawking temperature, as given in eqs.(|2j 
(129]) . depends on the values of parameters k and 70 , T = T(k, 70). Once 
these parameters are chosen, the coefficients a® and f\ of the metric functions 
expansion near the horizon can be calculated and then T determined from 
eq.([29]). 

According to the AdS/CFT correspondence p] -[3], properties of the dual 
3 dimensional field theory defined on the boundary can be read from the 
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Figure 3: Solution for the metric functions / and A as kHq grows from 
k = 0.1 (highest curve) to kHq = 1.2 (lowest curve) with 70 = 0.1 and 



Jo = 5.5. For this case n c 



1.24. 



asymptotic behavior of the solution of the system in the bulk. In this ap- 
proach the Hawking temperature T(k, 70) corresponds to the temperature of 
the d = 3 system on the boundary. 

Let us first consider the vacuum expectation value in the d — 3 field theory 
for the operator Ok, dual to the function K, associated with the magnetic 
field on the bulk. It follows from the identification (Ok) ~ K\ with K\ 
defined in eq.f lTH]) that K\ = K\(T) can be taken as an order parameter for 
the system in the border. As discussed for different models [I]- [9] one can 
interpret this result by stating that a condensate is formed above a black hole 
horizon because of a balance of gravitational and electrostatic forces. One 
can from eqs.( l2ip - (f2"2"j) see that the dimension A[0#] of the order operator 
is given by [H] 

3 1 

2 + 2 



A[0 



A" 



-^/l + 4r<L 2 



^+ 



(35) 



From the solution that we have found numerically we conclude that a 
finite temperature continuous symmetry breaking transition takes place so 
that the system condenses at a critical temperature T c , as can be seen from 
the behavior of Ki(T) for T w T c in figure 4. 

By fitting the curve we see that near T c one has a typical second order 
phase transition with power behavior of the form 

(0 K ) ~ K x oc (T - T c ) 1/2 as T ^ T c (36) 
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Figure 4: The K\ coefficient in the gauge field A{ asymptotic expansion as 
a function of temperature for J = 5.5 and kH = 0.2. To the right of the 
critical temperature (uncondensed phase) the solution corresponds to the 
trivial one, K\ = 0. The critical temperature for this choice of parameters is 
T c = 0.1798. 

The free energy T associated to the dyon solution can be calculated from 
the the on-shell Euclidean action action, 

J 7 = TSe\oti shell (37) 

We have numerically computed J 7 for the dyon solution and compared it with 
the free energy in the uncondensed phase which corresponds to the trivial 
solution ( 132]) . We plot in figure 5 the free energy density difference A J 7 , 
which can be seen to be continuous at T = T c . 

Other quantities of physical interest like the superconducting charge den- 
sity can be also computed following the same approach. The main conclusion 
we can draw form the calculations is that, concerning regions I and II in pa- 
rameter space, the results obtained by taking into account the back reaction 
of the matter on the geometry are similar to those obtained in ref. [9] where 
the Yang-Mills-Higgs system is studied in an AdS black hole background. 

On the other hand, the main difference with the fixed background case 
is that here we find a region in parameter space where a second horizon 
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Figure 5: The free energy difference between the condensed and the uncon- 
densed phases. The figure corresponds to kH = 0.2. To the left of the 
corresponding critical temperature J-dyon < -^trivial- 

arises and no non-trivial solution exists. This implies that one cannot arbi- 
trarily increase Newton's constant and the cosmological constant and have a 
nontrivial solution leading to a phase transition on the border. 

It is important to mention that in the gauge-gravity duality the extra 
dimension r does not describe an additional dimension but parameterizes 
the energy scale in the 2 + 1 dual field theory. So, the results described 
above indicate that the theory in the boundary corresponds to a system 
undergoing a superconductor phase transition as shown by the non-trivial 
order parameter that reflects the existence of a condensate. 

Let us recall at this point that, as explained before, the SU(2) gauge 
invariance of the Lagrangian in the bulk is spontaneously broken to a U{1) 
local symmetry due to the non-trivial condition one imposes to the Higgs 
field at infinity. Such gauge symmetry in the bulk corresponds to a U(l) 
global symmetry in the dual field theory, broken by the condensate which 
also breaks the rotation symmetry in the spatial plane. A combination of 
the two, the diagonal subgroup, is preserved. Concerning the fate of the 
corresponding Goldstone boson, in a superconductivity context one should 
start by coupling the theory on the boundary to a dynamical photon. It can 
be seen that the Goldstone boson is eaten by the photon, which becomes 
massive [TJ. 
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As it is well known it is the SU(2) — > U(l) symmetry breaking pattern in 
a gauge theory coupled to a Higgs field that allows the existence of regular 
monopoles and electrically and magnetically charged dyons. One can think 
that the further breaking of the local U(l) to the global U(l) in the boundary 
could be at the origin of the superfluidity/superconductivity phase regime in 
the presence of the dyons. Electric repulsion between dyons and the charged 
horizon compensates gravitational attraction and this is the reason why the 
superconducting layer floats above the horizon. 

Let us end this section by mentioning that it is possible to perform a 
simple analytic treatment leading to qualitative results in agreement with the 
detailed numerical analysis discussed above. The method starts by changing 
the radial coordinate to z = r^jr so that the horizon is located at z = 1 and 
the boundary at infinity at z — 0. One then proceeds to expand functions 
H, K and J near z = 1 and near z = and matching them at z — 1/2 by 
imposing continuity and smoothness. Using this method, originally proposed 
in this context in [27] for a simpler model, one can determine the coefficients 
of both expansions and compute relevant physical quantities. An analysis of 
this kind applied for a general family of models is in progress and we can 
advance that one has a good agreement in the case of the model discussed 
here. 



4 The Wilson loop 

As stated in the introduction, in low dimensional cases (d < 3), gauge fields in 
asymptotically AdS^+i spaces can induce on the boundary a field theory with 
dynamical gauge fields [UJ-[2D]. On this basis, following the gravity/gauge 
duality approach [2T]-j22], we shall compute the Wilson loop for such gauge 
fields in 2 + 1 dimensions in terms of the Nambu-Goto action for the string 
world sheet defined from the metric solution found in the precedent section. 
Since such solution corresponds to an asymptotically AdS black hole, our 
calculation corresponds to finite temperature |23]-|25j. 
We start from the metric ([H]) 

ds 2 = G tt dt 2 + G yy dy 2 + L 2 y 2 {d6 2 + sin 2 0d<p 2 ) (38) 



where 



G u = -A 2 (y)f(y) , G yy = J- (39) 

f(y) 
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and consider a string world-sheet parameterized in the static gauge according 
to 

t(r,a) = r, cp(r, a) = — ?-= , 9(r,a)=9, y(r,a) = y{cr) (40) 

Lsmu 

where 9 is a constant and y(a) defines the embedding, with a G (—1/2, 1/2). 
The quark q and anti-quark q are located at points (ip = a/2, 9 = 9) and 
(if = a/2, 9 = 9) respectively, separated a distance / along the "parallel" 
9 = 9. 

The resulting induced metric h a p(a) = Gmn(X) d a X M (a)dpX N (a) 
reads 

h = G tt (da ) 2 + (I 2 y(a) 2 + G yv y'(a) 2 ) (da) 2 (41) 

so that the Nambu-Goto action SV^ = T s J da a/— det h aj3 (a) can be writ- 
ten as 

i i 

/+2Z 
dx ^/F(u) 2 + G(u) 2 u'(x) 2 (42) 

with u(x) = y(Lx/l), F(u) = uA(u)fJ,(u)z and G(u) = A(u). 
A minimal action configuration has energy 



E(l) = T S L \F(u m ) l - + 2 jT" du G(u) Jl - ^^) (43) 

where u m is the minimum value taken by u and the distance I between the 
quark and the antiquark is given by 

1 = 2 T F{Um)G{U) ' (44) 

i„„ v/F( M ) 2 (f («) 2 - f (u, n y 



Energy E as given by eq. (l43|) diverges linearly for u^ — > oo. However 
by subtracting the bare mass m q of the two quarks, each one represented 
by long strings extended along the it-direction and puncturing the boundary 



Uu "'OO ^^ *k —I— or; 



2m q = 2T s L du G(u) (45) 
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we get a finite binding energy in the form 

E (l) = E(l) -2m q = T s L (F(u m ) l --2K\^ 



K [ l j)=[°° duG(u) U-Jl-^L. |+ / duG(u) (46) 



";, 



where Uh = 1 is the position of the horizon. We can re-write the quark- 
antiquark potential in the form, 

M = T, U„) - 2 ^f-\ (47) 

The quark- ant iquark binding energy E = E(l,T) can be calculated nu- 
merically by eliminating u m between eq.( l44|) and eq.( l4*6l) . Before doing this, 
let us analyze the behavior of each one of these functions in terms of u m . 

One can see in figure 6 that there exists a maximal distance l max between 
the quark-anti-quark pair, a phenomenon already found in the study of the 
finite temperature Wilson line in the large N limit of U(N) M = 4 Super 
Yang-Mills theory in 4 dimensions [23]- [25]. 

Concerning E = E(u m ), the generic behavior in a wide range of param- 
eters values is the one depicted in figure 7. Note that E = for a value u c m 
satisfying u c m > u r ^ LX leading to l c < l max . Thanks to these relations, valid 
in the whole parameter range that we have studied, the possibility of multi- 
valuation of the energy E as exposed by the figure is avoided: as I grows from 
/ = and reaches / c , the energy of the string configuration is the same as the 
energy of a pair of free quark-antiquark. From there on the answer given by 
formula (j47j) does not correspond anymore to the lowest energy which is in 
fact the free quark energy E = 0. 

Using eqs. f)44p - (14T|) one can determine the energy £ as a function of 
the distance I between the quark and the antiquark. The result is presented 
in figure 8 where we see that there exists a critical distance l c at which the 
quarks become free so that from there on one has a screening behavior. 

We have analyzed the behavior of l c as a function of the temperature 
finding that such distance shortens as the temperature grows, a behavior 




Figure 6: The distance I between the quark-antiquark pair as a function 
of u m , the minimum value of the variable u. Note that I has a maximum 
possible value, l max . The curve corresponds to the choice 70 = 0.1, k,Hq = 0.5 
and Jo = 5.5; for these parameters u™ x = 1.15 and l max = 0.90 (in units of 
eH ). 



E 




-0.8 1 



Figure 7: The energy of the quark-antiquark pair relative to the free quark 
configuration as a function of u m , the minimum value of the variable u (70 = 
0.1, kH = 0.5 and J = 5.5). For these parameters one has (in units of eH Q , 



ii-; 



1.44 and l c = 0.79. 
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E 




Figure 8: The energy of the quark- ant iquark pair as a function of I (70 = 
0.1, kH = 0.5 and J = 5.5). The graph corresponds to a temperature 
T = 0.018. 

that can be taken as a signal that the screening phenomenon is produced by 
the thermal bath. A representative graph is given in figure 9. 



5 Summary and discussion 

One of the purposes of the present work was to study the d = 3 + 1 di- 
mensional Einstein- Yang- Mills-Higgs model and find classical solutions to its 
equations of motion in asymptotically AdS space. As it happens for asymp- 
totically flat space [I2]-[13] we have shown that one can construct topo- 
logically stable self-gravitating dyons with a metric resulting from the back 
reaction on the space-time geometry solution that corresponds to a black hole 
in asymptotically AdS space. We find such solutions in a certain domain of 
parameter space which corresponds to region I in figure 1. There is also 
a region where only the trivial solution, for which the Higgs field takes its 
v.e.v. in the whole space, the gauge field is a pure gauge and the metric cor- 
responds to a Schwarzschild-AdS black hole (region II). As for the behavior 
of the Higgs and gauge fields spherically symmetric solutions, they are very 
similar as those we presented in ref. [9] for the case of a Schwarzschild-AdS 
black hole background (no back reaction), with the horizon radius inside the 
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Figure 9: The distance l c at which the quark and the anti-quark become free 
as a function of temperature. (70 = 0.3). 

dyon core and the Higgs field rapidly taking its vacuum expectation value. 
Finally there is a critical line which marks the appearance of a double hori- 
zon solution like the Reissner- Nordstrom metric. Above such line, no solution 
exists. 

Using the gauge-gravity correspondence, we were able to study properties 
of the dual d — 2 + 1 quantum field theory on the boundary of the 3 + 1 
asymptotical AdS space where our classical solutions were found. Because 
those solutions correspond to a black hole metric, they have an associated 
Hawking temperature T = T(7o, n) with parameters 70 and n related to 
the cosmological and Newton constants respectively. Then, varying those 
parameters in the bulk classical solutions we were able to study the finite 
temperature behavior of the vacuum expectation values of operators defined 
on the boundary, identifying an order parameter and showing that a second 
order phase transition takes place. 

Finally, because we have taken into account the back reaction of matter 
on geometry we were able to use the resulting black hole metric to calcu- 
late the Nambu-Goto action for our black hole solution related through the 
holographic correspondence to the Wilson loop at finite temperature. In this 
way we were able to evaluate the binding energy E = E(l,T) of two exter- 
nal charges in the boundary (a quark- ant iquark pair). We found a inverse 
power law behavior for E as a function of the distance between quarks, up to 
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a maximal distance l c where E vanishes indicating that the quarks become 
free in a typical screening process which can be attributed to the thermal 
bath. It should be remarked that this scenario is qualitatively the same al- 
ready found on the pioneering works on finite temperature Wilson loops in 
the large N limit of U(N) M = 4 supersymmetric Yang-Mills theory in 4 
dimensions |23j-[25j. 

Let us note at this point that, following the approach applied here, one can 
establish the existence of confinement for Yang-Mills theory in three dimen- 
sions by computing Wilson loops at finite temperature using the string/gauge 
correspondence (see ([28]) for a detailed calculation). We have found, in con- 
trast, a screening behavior so that we can discard the possibility that the 
model in the boundary could be related directly to 2 + 1 Yang-Mills. It is 
however well known that when massive fermions are added to a 2 + 1 gauge 
theory the effective action resulting from integration of fermions induces a 
Chern-Simons term and, as a result, confinement is destroyed (see for exam- 
ple [29]). Then, one can speculate that the theory on the boundary could 
include a Chern-Simons term. This idea may be tested by going beyond the 
classical approximation on the gravity side which would require to embed 
the gravity model into string theory. We hope to present a thorough study 
on these issues in the future. 
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